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Abstract 

The aim of this paper is to bring into the picture a new phenomenon in the theory 
of orthogonal matrix polynomials satisfying second order differential equations. The 
last few years have witnessed some examples of a (fixed) family of orthogonal matrix 
,S^ • polynomials whose elements are common eigenfunctions of several linearly independent 

j^ I second order differential operators. We show that the dual situation is also possible: 

there are examples of different families of matrix polynomials, each family orthogonal 
with respect to a different weight matrix, whose elements are eigenfunctions of a common 
second order differential operator. 

Theseexamples are constructed by adding a discrete mass at certain point to a weight 
matrix: W = W + StoM{to). Our method consists in showing how to choose the discrete 
i^Vj I mass Af(io), the point Iq where the mass lives and the weight matrix W so that the 

C*^ ' new weight matrix W inherits some of the symmetric second order differential operators 

associated with W. It is well known that this situation is not possible for the classical 

scalar families of Hermite, Laguerre and Jacobi. 

T-H ■ 

"^ ■ 1 Introduction 



The theory of matrix valued orthogonal polynomials starts with two papers by M. G. Krein in 
^ ■ 1949, see [Kl, K2]. A sequence of orthonormal matrix polynomials {Pn)n can be characterized 

^ as solutions of the difference equation 

tPnit) = An+lPn+lit) + B„P„(t) + A*^Pn^iit), n = 0, 1, . . . , 

where A^ and B^ are, respectively, N x N nonsingular and Hermitian matrices, and initial 
conditions P_i = and Pq nonsingular. Each family {Pn)n goes along with a weight matrix 
W and satisfies / PndWP^ = Sn,ml- 
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Nevertheless, more than 50 years have been necessary to produce the first examples of 
orthogonal matrix polynomials {Pn)n satisfying second order differential equations of the form 

(1.1) P;^{t)F2it) + P;^it)Fiit) + Pnit)Fo=TnPnit), n = 0, 1, . . . . 

Here F2, Fi and Fq are matrix polynomials (which do not depend on n) of degrees less 
than or equal to 2, 1 and 0, respectively (see [DGl, GPT2]). Two main methods have been 
developed in the last five years to produce such examples: solving an appropriate set of 
differential equations (see [D2, DGl, DG4, Ddl]) or coming from the study of matrix valued 
spherical functions (see [GPTl, GPT3, PT]). These families of orthogonal matrix polynomials 
are among those that are likely to play in the case of matrix orthogonality the role of the 
classical families of Hermite, Laguerre and Jacobi in the case of scalar orthogonality. The 
complexity of the matrix world has however led to an embarrassment of riches if we compare 
with these very few scalar examples (see, for instance, the papers cited above). 

We point out here that the second order differential equation (1.1) for the orthonormal 
polynomials {Pn)n is equivalent to the symmetry of the second order differential operator 

(1.2) D = d^F2{t) + d^Fi{t)+d^Fo, d=^. 

at 

with respect to the weight matrix W, as far as the eigenvalues r„ are Hermitian (see [Dl]). 
The symmetry of D with respect to W is defined by J{PD)dWQ* = J PdW{QD)*, for any 
matrix polynomials P and Q. Here and in the rest of the paper we will follow the notation 
in [GT] for right-hand side differential operators. In particular, if P is a matrix polynomial 
and D a differential operator as (1.2), by PD we mean 

PD = P"{t)F2{t) + P'{t)Fi{t) + P(t)Fo. 

As more families of orthogonal matrix polynomials satisfying second order differential 
equations become available many new and certainly interesting phenomena are being dis- 
covered. These phenomena are absent in the well known scalar theory. One of these 
phenomena is that the elements of a family of orthogonal polynomials {Pn)n can be com- 
mon eigenfunctions of several linearly independent second order differential operators. See 
[CG, D2, Ddl, DL, Gdl] for some examples illustrating this phenomenon. 

The purpose of this paper is to show what one can call the dual situation to that de- 
scribed in the previous paragraph. It turns out that there also exist second order differential 
operators as (1.2) having the elements of infinitely many different families of orthogonal ma- 
trix polynomials as eigenfunctions. More precisely, we will show examples of a fixed second 
order differential operator D as (1.2) for which there exist infinitely many weight matrices 
1^75 7 > 0, such that their corresponding monic orthogonal matrix polynomials (i-*„ 
eigenfunctions of D (with equal eigenvalues) 

Pn,^D = TnPn,f, n = 0, 1, . . . , 7 > 0, 

where D and r„ do not depend on 7. 



We give a simple but fruitful method to find such annoying examples (Section 2) and 
show a collection of instructive examples (Section 3). 

Our method itself is again a surprise if one compares with the situation in the scalar 
case. We first take a weight matrix W having several linearly independent symmetric second 
order differential operators. And then we add to W^ a Dirac distribution 5(gM(to), where the 
real number to and the mass M(to) (a Hermitian positive semidefinite matrix) are carefully 
chosen. We show in Section 2 that, fixed Iq, we can produce, under certain mild hypothesis, a 
second order differential operator D symmetric with respect to W and a positive semidefinite 
matrix M(fo) such that D is also symmetric with respect to any weight matrix of the form 
W + 7(5(0 M(to), 7 > 0. We illustrate with some examples, in Section 3, that the choice of 
the point to is not related to the support of W but to certain properties of the differential 
coefficients F2{t), Fi{t) and Fq evaluated at to. 

The situation is absolutely different to that of the scalar case. When we add a mass point 
to any of the classical weights of Hermite, Laguerre and Jacobi, the existence of a symmetric 
second order differential operator automatically disappears. Only when to is taken at the 
endpoints of the support, one eventually gets the symmetry of a fourth (or even larger) order 
differential operator which, of course, is not symmetric with respect to the original weight. 
This arises for the particular cases of the Laguerre weight e~* in (0, +oo), for the Legendre 
weight 1 in (—1, 1) and for the special case of the Jacobi weight (1 — t)" in (0, 1), raising the 
so called Laguerre type weight e~* + M6o, Legendre type weight 1 + M((5_i + 6i) and Jacobi 
type weight (1 — t)" + M6o, respectively (see [LK]). 

For the benefit of the reader, we conclude this Introduction displaying one of our examples. 
Consider the weight matrix 

W^a(t) = e-*' (^ +j'*' "/V teM, aeR\{0}. 

The linear space of differential operators of order at most two having the orthogonal poly- 
nomials with respect to Wa as eigenfunctions has dimension five. A basis is formed by the 
identity and four linearly independent operators of order two (see Section 6 of [CG]). We 

show in Section 3.1 that the weight matrices Wa^-y = Wa + jSq ( ^ i ) , 7 > 0, share the 

following symmetric second order differential operator: 



1 1 



D.^aM':,"' ";tI?V'''(~'T" 2a + 2(2 + a=)«,_^g„ 



^-1 2^' 
, " 1 



As we pointed out above, this means that each orthogonal monic family {Pn,a,'y)n, 7^0, 
with respect to Wa,'y satisfies the same second order differential equation, namely 



^n,a,'y^a — -L n,a^n,a,^^ ^ — U, i, . . . 



where 



2n + 1 



rn,a = I ' 4 ' «' 

^ -2n + l 

Notice that neither Da nor r„ „ depend on 7. 

We wiU also see in Section 3.1 that the choice of the point to where the discrete mass is 
added can be located at any real number. There always exists a Hermitian positive semidef- 
inite matrix M{to) such that the weight matrices Wa + 'jSt^M^to), 7 > 0, share a common 
symmetric second order differential operator (see (3.2)). 

2 The main result 

This section is devoted to present a set of constraints to guarantee the symmetry of a dif- 
ferential operator of any order with respect to a weight matrix modified by adding a Dirac 
distribution at an arbitrary point. Before that, we need some definitions and previous results. 
Throughout this paper W will denote an N x N weight matrix supported in the real line. 
Then, we can consider for any pair of matrix polynomials P and Q the Hermitian sesquilinear 
form 



{P,Q) = j P{t)dW{t)Q*{t), 



where Q*{t) denotes the conjugate transpose of Q{t). 

As far as W is non-degenerate, we can produce a sequence of orthogonal matrix polynomi- 
als {Pn)n with deg Pn = n and nonsingular leading coefficient such that {Pn, Pm) = ^n^n,mi 
with A„ nonsingular. If A„ = /, we say that the sequence {Pn)n is orthonormal. 

The non-degenerateness of W can be defined in several ways. For instance, by assuming 
that J P{t) dW{t) P*{t) is nonsingular if the leading coefficient of the matrix polynomial P 
is nonsingular. A simple way to check that W is non-degenerate is the following: consider 
the absolutely continuous part of W with respect to the Lebesgue measure and write W{t) 
for its Radon-Nikodym derivative. Then W is non-degenerate if W{t) is positive definite in 
an interval of the real line (for a much more complete introduction to matrix orthogonality, 
see [DG3] and references therein). 

If one is considering possible applications of orthogonal matrix polynomials, it is natural 
to concentrate on those cases where some additional property holds. For instance, in [Dl] 
one of the authors raised the problem of characterizing weight matrices whose orthonormal 
matrix polynomials are common eigenfunctions of some symmetric right-hand side second 
order differential operator D as (1.2) with Hermitian eigenvalues. We say that a differential 
operator D is symmetric with respect to a weight matrix W if 

{PD,Q) = {P,QD) 

for any pair of matrix polynomials P and Q. We recall that we are following the notation 
in [GT] for right-hand side differential operators. We already mentioned in the Introduction 



that in the last few years a large class of families of weight matrices W has been found having 
symmetric second order differential operators as (1.2). 

The condition of symmetry for the pair made up of a weight matrix W and a differential 
operator D/. of order k can be established in terms of a set of difference and differential 
equations relating W and the coefficients of D^. Indeed, if we write the right-hand side 
differential operator D^ of order k as 

k 

(2.1) Dk = Y,d'F,{t), d = j^, 

i=0 
where Fi{t),i = 0, . . . ,k, are matrix polynomials of degree less than or equal to i, 



j=0 



^NxN 



and denote by fin, n = 0, 1, . . ., the moments of the weight matrix W, i.e., /in = J t'^''dW{t), 
then we have the following 

Theorem 2.1. For a weight matrix W the following two conditions are equivalent: 

1. The operator D^ is symmetric with respect to W . 

2. For n > I, the following k + 1 sets of moments equations hold 



j=0 

where 



Bl = Y,FLlin-^, l = 0,...,k. 



i=0 



Moreover, suppose the weight matrix W = W{t)dt has a smooth density W{t) with respect to 
the Lebesgue measure which satisfies the boundary conditions that 

(2.3) E(-l)'"^'''"'(^70(^^--^)^'"'"^' ^=1'---'^' l = 0,...,k-p, 



1=0 



should have vanishing limits at each of the endpoints of the support of W , and the following 
k + 1 matrix differential equations hold 



(2.4) ^(_i)'=-(^ ']{Fu^,.wf-'^'^=W.Ft, l = 0,...,k. 



j=0 



Then the differential operator D^ (defined in (2.1)) is symmetric with respect to the weight 
matrix W. 



(In (2.2) we are using the well known notation for the falling or bounded factorial (x)„ 
defined by (x)„ = x{x — 1) • • • (x — n + 1) for n > 0, (a:)o = 1.) 

Proof: The first and second part are shown to be equivalent in Proposition 4 in [Ddl]. 
Likewise, last part can be found in Theorem 5 in [Ddl], using integration by parts (taking 
into account the boundary conditions (2.3)). 

D 

In particular, the symmetry equations (2.4) for differential operators of order two appeared 
for the first time in [Dl] when FqW = WFq and later in [DGl] and [GPT2]. As we remarked 
in the Introduction, most of the examples which appeared in the last years have been obtained 
solving these equations, while some others came from group representation theory. Moments 
equations (2.2) turn out to become a suitable tool in all examples throughout this article. 

We now present the result that we will use to generate examples of second order differential 
operators having infinitely many families of orthogonal matrix polynomials as eigenfunctions. 
The idea is to find certain constraints which guarantee the symmetry of a differential operator 
with respect to both W and a new weight matrix obtained from W by adding of a Dirac 
distribution at one point. 

Let H^ be a weight matrix and consider the following weight matrix 

(2.5) W{t) = W{t) + 6t,{t)M{to), 

where (5t,j(t) = 5{t — to) is the Dirac delta distribution or the "impulse symbol" introduced 
by P. A. M. Dirac in [D], which we will consider as a measure. The Hermitian positive 
semidefinite matrix M(to) depends on the point where the Dirac distribution is added. 

Weight matrices of the form (2.5) were considered in [YMPl, YMP2] (to study asymptotic 
properties of the corresponding modified Jacobi matrix) for a weight matrix W in the Nevai 
class, i.e., with convergent recurrence coefficients. 

The moments of W are related with the moments of W by the formula 

Jln= [ t'^dWit) =fln+ [ t^'K {t)M{to)dt = fin + t^M{to), n = 0, 1 



5 ^5 



Observe that in the special case of to = the only modified moment is the first one JIq = 
/xo + -^(0), and then pt„ = ^„ for n = 1, 2, . . .. 

The following theorem gives conditions for the symmetry of a differential operator Dk 
with respect to the weight matrices W and W + St(^M{to). 

Theorem 2.2. Let D^ he a differential operator of order k as in (2.1). Let W he a weight 
matrix. Assume that associated with the real point to € IR there exists a Hermitian positive 
semidefinite matrix M(to) satisfying 

(2.6) F,(to)M(to) = 0, j = l,...,k, 

FoM{to) = M{to)F^. 

Then the operator D^ is symmetric with respect to W if and only if it is symmetric with 
respect toW = W + 6t^,M{to). 



Proof: Recalling definitions around (2.2) ior W = W + 6tQM{t()), we produce 
Bl = Y.^LJ^r.-i = Bl + f^-'Fi{to)M{to), l = 0,...,k. 

i=0 
Using conditions (2.6) for j = 1, . . . ,k, we obtain 

5o = SO+tJFoM(to), Bl = Bl, l = l,...,k. 

Consequently, this shows that equations (2.2), I = 1, . . . ,k, are just the same for W and W. 
For / = 0, equations (2.2) for W and W are, respectively: 



fc-i 

j;(n),_.i?r^ + i?o = (i?0)*, 

i=0 
k-l 

Y;^{n)k^,Bt' + t^FoM{to) = {B^J* + toM{to)F^. 

i=0 

The last condition in (2.6) shows again that those equations are the same for W and W. 

D 

We emphasize that Theorem 2.2 for A^ = 1, i.e., for scalar orthogonality, implies that 
either M = or there exists a common zero for all coefficients of the differential operator. 
For instance, for k = 2, the existence of that common zero is not possible for the classical 
families of Hermite, Laguerre and Jacobi. 

Let us notice that once we generate W, Di^, to and M(to) satisfying constraints (2.6) 
of the theorem above, we can produce not only one single weight matrix for which D^ is 
symmetric, but also infinitely many different weight matrices for which D^ is symmetric as 
weh. If Theorem 2.2 fulfills for W+6toM{to) then automatically also fulfills for W+-f6toM{to) 
where 7 is any positive real number. 

For a weight matrix W, the constraints (2.6) mean that by adding a Dirac distribution to 
W, the chance that W and W + dt^M^to) share a symmetric differential operator of order k 
increases with the number of linearly independent symmetric differential operators of order 
k for W. In fact, all the examples we show in the next section are built from a weight 
matrix W having several linearly independent second order differential operators. Sometimes 
we can add the Dirac distribution at any real number: in our examples, this happens when 
the number of linearly independent symmetric second order differential operators is four. 
Sometimes we lose this freedom and we have to add the Dirac distribution at an specific real 
number: in our examples this happens when the number of linearly independent symmetric 
second order differential operators is just two. In this case, we have to locate the mass at one 
of the endpoints of the support. Both kind of examples reveal that the choice of the point to 
where the Dirac distribution is located depends more on the matrices i^2(io)) ^1(^0 ) ^ind Fq 
than on the support of W. 



3 Examples 

In this section we exhibit a collection of instructive examples. The first three examples of 
(2 X 2) weight matrices W (supported in (— oo,+oo), (0, +oo) and (0,1), respectively) have 
the property that they provide four linearly independent second order differential operators 
having a fixed family of orthogonal matrix polynomials with respect to W as eigenfunctions. 
As we mentioned above, for any real number Iq we can find a positive semidefinite matrix 
M(io) and a symmetric second order differential operator D as in (1.2) satisfying constraints 
(2.6). According to Theorem 2.2, the operator D will be also symmetric with respect to 
W^ + 7'^to^(*o),7>0. 

The last two examples deal with {N x N) weight matrices W (supported in (0, +oo) and 
(0, 1), respectively) with at least two linearly independent symmetric second order differential 
operators. This forces us to locate to at one of the endpoints of the support of W with a 
mass M carefully chosen. 

More examples will appear in the PhD dissertation of one of the authors [dl] . 

3.1 Wait) = e-*'e^*e^** with ^ = (q qY a G M \ {0} 

This weight matrix was introduced for the first time in Section 5.1 in [DGl] (for the general 
size NxN). It was deeply explored in [DG2] and a set of generators of second order differential 
operators can be found in Section 6 in [CG]. 
By expanding the exponential, we find that 

(3.1) VF„(t) = e-*'e^*e^** = e-*'('^+j'*' "^^V t G M, aGM\{0}. 

Using (2.4) (for /c = 2), we find an expression for the 5-dimensional (real) linear space 
of symmetric differential operators of order at most two with respect to Wa- Then, for a 
fixed real number to, we solve the equations (2.6). In this case, we find that the following 
differential operator 

(3.2) Da,to = d^F2{t) + d'F,{t) + 5°Fo, 

where 

P.,^ f-Clto+ato-at -l-{aHo)t + aH^ 

p (,-, _ ^-2« + 2e^t„i -2to - 2<^,^, + 2(2 + a^)t\ 
'^' V 2to 2{C,t,-ato)t J 

/ 3- to 2 + 0^ 

- -e^ -2-° 



and the Hermitian positive semidefinite matrix M(to) 



M{to) = M{a,to) = (^^'r^ ^Y 



where 






satisfy the constraints (2.6). 

This differential operator can be obtain by a hnear combination of the second order 
differential operators introduced in [CG] (Dj, i = 1,2,3,4), namely 

o.... = (-&. + ^)^-S,„A ->. + >. 

It is easy to verify, using ^a,to^a,tQ + 1 = 0, that the coefficients of -Da,to evaluated at to 
satisfy the conditions (2.6). Since Da, to ^^^ been chosen to be symmetric with respect to the 
weight matrix (3.1), Theorem 2.2 implies that Da^to is also symmetric with respect to any of 
the following weight matrices: 

Wa,torf{t) = Wa{t) + j6t,{t)M{a,to), j>0. 

As a consequence, the sets of monic orthogonal matrix polynomials iPn,a,to,-y)n with re- 
spect to Wa,to,7' 7 — 0' ^^^ eigenfunctions of a common second order differential operator, 
namely -Da, to; defined in (3.2). That is, 

^n,afy,'y^afy — -L n,a,to^n,a,ti^,^^ ^ — 'Jj -L, . . . , 

where 

1 n,a,to - 4 

V ^ eJ,„(2n-l)-2nato- 

Notice that neither -Da,to ^lo^ '^n,a,to depend on 7. 

We would like to make an special mention to the weight matrix W{t) = e~* e e"^** where 
A is the 3x3 matrix 

/O a 0\ 
A= b\ , a,6ElR\{0}. 
\0 0/ 

Assuming that 

(3.3) a^h'^ = A{}?-a^), 

the dimension of the linear space of differential operators of order at most two having the 
orthogonal polynomials with respect to W as eigenfunctions is just three (if (3.3) does not 
fulfill, the dimension is two, while for the weight matrix Wa in (3.1) is five). A basis is formed 
by the identity and two symmetric second order linear operators (see [D2]). Our method, 
however, does not provide any common symmetric second order operator for W and any 
modification of W by adding a Dirac distribution. 



3.2 WaAt) = f^e-H^t^* with B=[\ JJ j , a G M \ {0} 

This weight matrix was introduced for the first time in Section 6.2 in [DGl] (for the general 
size A^ X N) and it was extensively studied in [DL]. Unlike the first example, its algebra of 
differential operators has not been studied in depth, but as we mentioned at the beginning 
of this section, it turns out that there are four linearly independent symmetric second order 
differential operators and the behavior of the algebra is quite similar to the one introduced 
in Section 3.1. 

By computing t^ for B = [ n ) ' ^^ ^^^ that 

.-UB,B' _ ,c.^.-t (t^ + aHt - 1)' a{t - r 



(3.4) VF,,«(t)=t°e-'t^t^ =re-M \^^^_^.' 1 j' *^(0'+°°)' « > "l- 

Proceeding as in Section 3.1, for a fixed real number to, we find that the following differ- 
ential operator (symmetric with respect to Wa a) 



where 



Da,c.,to = d^F2{t) + d^Fiit) + 9%, 



ato «'^o^^JC^- '+^"+y+°^ -(a + to + l)(l + a2^ 



to -ato) ' ~ \ C^ + a 

+ t 



o/O (l + a2)(l + a; 







-a(3to - 2 + (q + 1)2) \ /2aC^ - (1 + a^)(a^ + 2to + 3a; 
+ (a + 3)(^±_*o^)^ V +ato(a2-l)-(to + a + 3) 
to — a — \ C^{a + 1) — aato 

C± + a{to - 1) + ^fl^ (1 + a^){a^ + 3a - t^a + 2) + 2(q - V 
-C^ + aa 



l+g C^ «(fo-l) (a+2)(fo+a+l) , a(l+a) 



T+c? 2 2 2a + 1+a^ 



and the Hermitian positive semidefinite matrix M{tQ) 

-± \2 A± 



\ ''a, a, to / 



where 



±^ ± _ 1 (g^ + l)(to + a)-a^ + l± ^(g^ + l)(a2(to - a - 1)^ + (tp + a + Ip) 
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satisfy the constraints (2.6). 

Thus, Theorem 2.2 imphes that Da,a,to is symmetric with respect to any of the following 
weight matrices: 



W^a,a,to,7(*) = WaAt)X{id,+oo){t) + lK{t)M{a,a,tQ), t G 



a > -1, 7 > 0. 



Note that the discrete mass on the Delta distribution can be located in or out of the support 
of the original weight matrix (3.4). 



3.3 An example supported in (0, 1) 

The following weight matrix is a modification of the one introduced in [PT] (see also (3.8) 
below in this paper): 



(3.5) Wa,p,kit) = r(l - t)f^ 



kf+f3-k + l {p-k + l){l-t) 
(/3-fe + l)(l-t) {p-k + l){l-tf 



tG(0,l) 



where a, (3 > —1 and < k < (3 + 1. 

The difference between Wa^f^^k and the weight matrix introduced in [PT] (A^ = 2), is that 
Wa^j3^k enjoys four linearly independent symmetric second order differential operators while 
that in [PT] enjoys only two. The example in [PT] is related to a group-theoretical situation 
but this is not the case of our Wa,i3,k, as far as we know. 

Like in the previous examples, for a fixed real number to (except for to = —■^{c(+l3 — k + 2) 
and to = 1), no matter if it is located in or out of the support of the weight matrix Wafi,k-, 
we find a symmetric second order differential operator Da,/3,k,to with respect to Wa,p,k and 
a Hermitian positive semidefinite matrix M(to) satisfying the constrains (2.6). Hence, this 
operator Da,f3,k,to is also symmetric for any of the weight matrices Wa,/3,k + ^^to^ (to) j 7^0. 
Since the formulas for arbitrary a, (3 and k are unpleasant, we show here only one concrete 
example: a = 0, /3 = and k = 1/2: 



Fo 



il-t){-to + ti2to-3)) + 
-toil-tf 



t(l-t)(l-to) 



V9- 



2t + to — 2tot — tot 

t(l-t)(l-to) 



(t-l)(t-to) + 



/l2t-10 + 8to(l-t) + ^^5 ^^y ^^ g_t-Ato{t + 2) + 



Lp- 



\ 



(l-4to)(t-l) 



4(t - to) + 




(to-l)(4t-l 



/ (l-to)(l + 2y^) 

2(/9± 

1-to 



-to - 3 + 



1-to \ 



V 



-to + 



2ip^ 

:i-to)(i+2(^±) 



2lp^ 



2ip- 



,± 



and 



V? 



^(*o) - \^^± (^±)2 ; , 
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where 



± 2 - to ± V2io - 2to + 1 
•^ (to) = 



io + 3 

We need to impose to 7^ —3 to avoid singularities in M(to) and to 7^ 1 in f^ (because 
99^(1) = 0) to avoid singularities in the differential coefficients F2,Fi and Fq. 

3.4 An example of general size 

We consider here the weight matrix defined by 

(3.6) W{t) = t°e-*e^*ti-^t^-^*e^**, t G (0, 00), a> -1, 
where 

(3.7) J = Y,iN-i)Eu, A=Y,^iE^,i+l, ly^€^\{0}. i = l,...,N-l. 

i=l i=l 

Here we are using Eij to denote the matrix with entry {i,j) equal to 1 and otherwise. 

This weight matrix was introduced for the first time in [Ddl]. This example enjoys the 
special property of having symmetric odd order differential operators, a phenomenon that is 
not possible in the classical scalar theory. For more details, the reader should consult [Ddl]. 

It is proved in [Ddl] that (3.6) has always a symmetric second order differential operator 
given by 

Di = dhl + d^[{a + 1)1 + J + t{A - I)] + d^[{J + aI)A - J], 

where A and J are defined in (3.7). Note that there are A^ — 1 free parameters in Di. 
Assuming the following conditions on the parameters i^i, . . . , i'n-2'- 

i{N - i)4_i = (iV - l)uf + {N-i- l)vjvl_^, i = 1, . . . , AT - 2, 

the weight matrix W has another symmetric second order differential operator 

D2 = d'^F2{t) + d^Fiit) + d^Fo, 

where the coefficients are given by 

F2=t{J-At), 

Fi = ((1 + a)I + J)J + Y- t{J + {a + 2)A + Y* -AY + YA), 

N -1 

[J-iaI+J)A], 



4-1 



A and J are defined in (3.7), and Y = \ -E^i+i 1. Note that now the only free 



parameter is vn-\- 



Yl 



Let us call M = M(q;, z^i, . . . , i/jv-i) the Hermitian positive semidefinite matrix (and 
singular) 

AT ,max{i,j}-l , ,7V-max{Jj} , M\2 

«=e( n "'^"m ,^% n i^^-f^U«. 

«j'=l fc=min{ij} A;=l 

Let Z? be the following differential operator 

D = -{N-l)Di + D2. 
Then, evaluating the differential coefficients Fi{t), i = 0,1,2, oi D at t = 0, we observe that 

i^2(o) = o, 

TV 

Fi(0) = Y- J2{i -l){a + N-i + l)Eu, 

i=l 
(iV-l)(l + ^.2^_l), 



iy2 



-[J - {al + J)A]. 



Writing down the expressions and considering the bidiagonal structure of -Fi(O) and Fq, it is 
easy to conclude that the conditions (2.6) are satisfied. Also, the second order differential 
operator D is a linear combination of symmetric differential operators with respect to the 
weight matrix (3.6), so it is symmetric. Hence, Theorem 2.2 implies that D is symmetric 
with respect to any of the following weight matrices: 

W^it) = Wit) + jS{t)M, t G (0, +CX)), a > -1, 7 > 0. 

We illustrate the case of size 2x2 (i^i = a). The weight matrix Wy is 

„ iA(l + o2t) at\ fa\a + lf a{a + 1)\ 

Wy(t) = t^ [ at ij+^U(" + i) 1 J ^^' 

while the second order differential operator D is 

2,0 -at'\ it _(1±^H^±3)1\ Ja^ + l (l + a2)(a + l) 




3.5 The one step example 

This weight matrix was introduced for the first time in [PT]. An associated second order 
differential operator comes from the Casimir operator acting on matrix valued functions on 
the group SU(n + 1). The reader is advised to consult, apart from [PT], either [GPT3] or 
even better [GPTl] for a full account. 
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The weight matrix, for t £ (0,1), a, f3 > —1 and < k < f3 + l has the fohowing structure: 

(3.8) WaMt)=t"i^-tfZ{t), 

where 



N / N 



~'^^-^^E,,. 



Likewise, as in Section 3.4, we have an expression of two symmetric second order differential 

operators. A first one Di is given in [PT] (Theorem 3.4) and a second one D2 is given in [PR] 

(Theorem 3.1). Notice that in [PT, PR], left-hand side differential operators are considered. 

Let us call M = M{a, [3, k) the Hermitian positive semidefinite matrix (and singular) 

>^ /iV-l\ (N-l\ {a + l3~k + i + 1)(^-^) {a + l3-k + j + 1)(^-J) 



where the number (x)^"' denotes the rising or upper factorial defined by {xY"'' = x{x + 
l)---(x + n-l) for n > 0, (a;)(°) = 1. 

Let D be the following differential operator 

D = -{a - N + l)Di + D2. 

Thus, evaluating the coefficients Fi[t), i = 0,1,2, oi D at t = we observe that 

F2(0) = 0, 

N-l N 

Fi (0) = Y^ 3i(a + (3-k + i + l)Ei^i+i -^^{3 - k + i - 1){N - i + l)Ei^i-i 

1=1 i=2 

N 



+ Y^ 3[(i -l){a + (i-k + i) + {i-A){l3-k + i)]Eii, 
1=1 

N N 

Fo = -3(iV + A: - 1) J^(i -l){a + l3-k + i)Eii + ^(iV - i + 1){I3 - k + i - l)Ei, 



i-l- 
1=1 i=2 

Writing down the expressions and considering the tridiagonal structure of -Fi(O) and Fq, it 
is easy to conclude that the conditions (2.6) are satisfied. Also, the second order differential 
operator D is a linear combination of symmetric differential operators with respect to the 
weight matrix (3.8), so it is symmetric. Then, Theorem 2.2 implies that D is symmetric with 
respect to any of the following weight matrices: 

Wa,/3,k,^(.t) = WaMt)+l^(.t)M{a,^,k), tG(0,l), a,/3>-l, 0<fc</3+l, 7 > 0. 
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Like before, we illustrate the case of size 2x2. The weight matrix W^ = Wa,/3,k,'y is 



Pi kt + P-k + 1 (1 




[l-t){P-k + l) (1- 

while the second order differential operator D is 

2^0 \ f -p + k-1 (p-k + l){t-l) 

\t t{l-t)J^ \a + P-k + 2 {a + p + 4){l-t)-{k + 2] 

.ao/'O {l + k){P-k + l) 
^ \0 {l + k){a + P-k + 2] 
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